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ABSTRACT 

In addition to its intrinsic practical importance, 
nonlinear time delayed feedback control applied to 
lifting surfaces can result in interesting aeroelastic 
behaviors. In this paper, nonlinear aeroelastic response 
to external time-dependent loads and stability boundary 
for actively controlled lifting surfaces, in an 
incompressible flow field, are considered. The structural 
model and the unsteady aerodynamics are considered 
linear. The implications of the presence of time delays 
in the linear/nonlinear feedback control and of 
geometrical parameters on the aeroelasticity of lifting 
surfaces are analyzed and conclusions on their 
implications are highlighted. 

NOMENCLATURE 

a Dimensionless elastic axis position measured 
from the mid-chord, positive aft 
b Semi-chord 

c , k Damping and stiffness parameters of 1-DOF 
plunging airfoil, respectively 
c h ,c a ,Cn Damping parameters in plunging, pitching 
and flapping, respectively 
C(,v) Theodorsen's function 
C La Lift-curve slope, 2 n 
e Dimensionless leading edge flap position 
measured from the mid-chord, positive aft 
F a , F b Aerodynamic and time-dependent load vectors 
F c Control force 

g p , g v , g nc Proportional, velocity and nonlinear 

feedback control gains matrix 
g P , g v , g nc i > 9v > 9« c Proportional, velocity and 

nonlinear feedback control gains of 1-DOF 
airfoil, and their dimensionless counterparts, 
respectively 

G Control input matrix 
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/?,«,/? Plunging, pitching and flap displacements, 
respectively 

I a , Ip Mass moment of inertia per unit length of the 

wing-flap and of the flap about the elastic axis, 
and about the flap axis of rotation, respectively 
k h ,k a ,kp Stiffness parameters in plunging; 

torsional stiffnesses of the wing and flap about 
the elastic axis, and about the flap axis of 
rotation, respectively 
k Reduced frequency, cob/U x 
L a ,L h , L c Aerodynamic lift, time-dependent external 
load, active feedback control, respectively 
m , p Mass of the wing per unit of length and mass 

ratio, (= m/lp^b 2 ^, respectively 
M , K , B Structural matrices, 

M a , K a , B a Aerodynamic matrices 

r a , rp Dimensionless radii of gyration of the wing-flap, 

[l a / mb 2 ) , and of the flap, [i p j mb 2 ) , 
respectively 

s Laplace transform variable 

S a , x a Static unbalance about the elastic axis and its 
dimensionless counterpart, S a I mb 
S p , x b Static unbalance about the flap axis of rotation 
and its dimensionless counterpart, S p /mb 

t , cr ; r Time and dummy time variables, respectively; 

dimensionless time variable, tU x /b 
Tj Theodorsen’s constants 

U K ,LFreestream speed and its dimensionless 
counterpart V =U X /co 0 b , respectively, 
x Plunging, pitching and flap displacement vector 

£t,£ a ,£p Structural damping ratios in plunging 

(= c h / 2mco h ) > pitching (= c a /2l a co a ) , and 
flapping {= c p 121 p(o p), respectively 
4" , a Plunging and pitching displacement quantities 
p x Air density 

t t , Tj Time delay and dimensionless time delays, 
tjU x /b ,i = 1,4 , respectively 
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<j>[r) Wagner’s function 

a> h ,co a ,co a Uncoupled frequencies in plunging, 

(k h /mf 12 , pitching, (k a /I a ) 1/2 ,and flapping 
(kp jl 12 , respectively 

INTRODUCTION 

ast and recent literature on aeroelasticity is mostly 
devoted to linear models and to harmonic solutions. 
Often experimental results are interpreted by assuming a 
linear behavior of the physical model. Recently, special 
emphasis has been placed on the role of nonlinearities 
on aeroelastic instabilities, and for the most part, these 
studies have focused on the qualitative nonlinear 
behavior of open loop aeroelastic systems. Still, it is 
rather accurate to say that currently there is a lack of 
research work on closed loop dynamics of aeroelastic 
systems. The nonlinear aspects that we are addressing in 
this research are those which arise in the description of 
the feedback delayed control. 

This study can lead to new qualitative results in the 
areas of flutter instability boundary and aeroelastic 
response to gust and blast loadings. 

As a bottom line, it is imperative that the occurrence of 
flutter phenomena be suppressed, as to avoid the 
catastrophic failure of the structure [1]. These facts 
emphasize the importance of developing proper 
methodologies for the active control of structural 
systems, enabling one to raise the flutter speed, to 
enhance the aeroelastic response (attenuating excessive 
vibrations), and convert the unstable LCO, in which 
case the flutter boundary is catastrophic, into a stable 
LCO, for which case the flutter boundary is benign. 

The determination of the stability boundary of 
linear/nonlinear actively controlled aeroelastic systems, 
where the presence of the unavoidable time delays 
between controller and actuators is included, constitutes 
an important practical problem. In fact, the actuators 
may input energy at the exact moment when the 
controlled system does not need it [2-4]. These delays 
can be very detrimental in the sense of deteriorating the 
control performance and can even cause irregular 
motions, producing instability of the aeroelastic system. 
However, there are cases where those delays are used to 
control chaotic motions [5]. 

For a better understanding of the challenging problem 
related to the nonlinear delayed feedback control, the 
model of a 2-D wing section was considered and its 
effect on the aeroelastic response and flutter has been 
investigated. 

The methodology used in this work is based on Volterra 
series and indicial functions in conjunction with a 
feedback control [1]. Volterra’s functional series 
technique was proven to be an efficient tool in the 


solution of various nonlinear aeroelastic problems [4]. 
The Volterra series approach can also be used toward 
the formulation of the stability criteria for systems 
featuring time delays. 

BACKGROUND 

The problem of controlling unstable motion is an 
important subject in aeroelasticity. The renewed interest 
of aeroelasticians in this subject started, however, with 
the observation that a large number of unstable periodic 
orbits, (limit cycle oscillations LCO), embedded in 
chaotic attractors can be stabilized by weak external 
forces. 

Two main methods for controlling unstable motions 
have been established. The first one was developed by 
Ott, Grebogi, and Yorke [6]. This method, based on the 
invariant manifold structure of unstable orbits, is 
theoretically well understood, but since it is difficult to 
apply to fast experimental systems, it has a limited 
interest. Another approach is due to Pyragas [7]. This 
method uses time-delayed controlling forces; it can 
easily be applied to real experimental situations, but so 
far the control mechanism has been poorly understood 
from a theoretical point of view. Considerable research 
has been done for more than three decades on various 
aspects of dynamical systems with delayed factors in the 
state variables and/or control inputs [8-12]. 

Various stability criteria and numerical approaches have 
been presented in recent years, see [13-14] and 
references cited therein. Consequently, time-delayed 
feedback has been applied widely [3,15-17]. In Ref. [8] 
an investigation of the stability and chaos for wheel 
suspensions was presented. 

For aeroelastic systems in Ref. [5] the time-delayed 
feedback was applied toward the control of the chaotic 
motion of a 2-D lifting surface, with cubic pitching 
stiffness and linear viscous damping, using the feedback 
control method of Pyragas [7]. 

Stability analysis has been conducted in Ref. [2] for a 
linear, damped SDOF system with time delays in the 
displacement and in the velocity feedback. 

With the exception of [5], the use of time-delayed 
feedback in aeroelasticity has been very limited. This is 
due to the fact that the characteristic equation of the 
delayed system is transcendental, i.e. has infinite 
number of roots, so it is neither possible to solve for its 
roots, nor to easily find approximate solutions [2,12]. 
Moreover, the aeroservoelastic problem is extremely 
complex. For this reason, as a first step toward the 
nonlinear analysis, the stability of linear differential- 
difference aeroelastic equations has to be studied [3]. 

The present study can provide broad information and 
answer some basic questions, such as whether the 
aeroelastic stability is affected by the presence of delays 
that appear in the feedback, and whether the system 
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stability is robust with respect to small variations of the 
feedback gains. 

As reported in [4], multiple degree-of- freedom 
aeroelastic systems, including structural and 
aerodynamic nonlinearities, can be investigated via a 
combined Volterra series [18,19] and Indicial functions 
technique [1,4]. Originally, the methods of Volterra 
series and Volterra kernel identification were developed 
to identify the nonlinear behavior in electrical circuits 
[19]. In [1] Volterra’s series approach has been applied 
to the open/closed-loop aeroelasticity of airfoils. It was 
shown that the method provides opportunities for 
developing a unified and efficient way to address 
problems of nonlinear aeroelasticity. By performing an 
analytical linear stability analysis of the aeroelastic 
system via the use of the first order Volterra kernel, it is 
determined which class of orbits is accessible to time- 
delayed feedback control methods. Explicit expressions 
for important quantities like the critical time-delays and 
control gains or the dependence of the transient 
behavior on the control parameters, are derived. In this 
paper the stability boundary of a reduced order open/ 
closed-loop aeroelastic system incorporating a nonlinear 
time delayed feedback control is presented. The goals of 
designing such a system consist of controlling the 
aeroelastic response behavior, increasing the flutter 
speed and converting the catastrophic flutter boundary 
into a benign one. 

ANALYTICAL DEV ELOPMENTS 

The first step towards the modeling of an open/closed- 
loop aeroelastic system with nonlinear time delayed 
feedback control via the Volterra series approach is to 
determine the aeroelastic kernels. For the purpose of the 
present analysis, the approach presented in Refs. [1,4] 
has been modified. The basic assumptions and a 
detailed procedure are presented in Ref. [4]. For 
exhaustive treatments of the Volterra series concept 
applied in structural dynamics, the interested reader is 
referred to Ref. [20]. 

The aeroelastic kernels including control effects are 
derived in terms of the structural parameters, unsteady 
aerodynamics, proportional (PFC) and velocity (VFC) 
feedback control gains and feedback delays. Based on 
these kernels, the time histories and flutter boundary of 
the open/closed loop delayed aeroelastic system are 
obtained. The determination for each specific flight 
condition of the corresponding linear and nonlinear 
kernels of the Volterra series is required [4]. The 
open/closed-loop aeroelastic governing equation of an 
airfoil featuring plunging-pitching-flap deflection 
motion and subjected to external time -dependent loads 
can be expressed as 

M At) + + K,x(t ) = nr' [F a (t) + F b (*)] + Gu(t) (1) 


where x(t) = [h(t),a(t), p(t) Y , u(t) is the control input 
(for a 3-DOF, for example, a torque applied at the flap 
[21]). The unsteady aerodynamic loads are represented 
by 

F.(*) = M a x(/)+ B a x(/) + K.x(f ) + F c (?) (2) 

The significance of the other parameters is well known, 
see [22,23]. As a remark, a closed-loop system can be 
seen as an open-loop system where the transfer function 
includes the feedback control. In Eq. (1) the state 
feedback control with delay can be expressed in the 
form: 

Gu(f -t)= g p x(t - r)+ g v x(t - r) + g nc x 3 (f - r) (3) 

where g p ,g v ,g nc are the feedback gain matrices for the 

displacement, velocity, and the nonlinear term, 
respectively. Since the aeroelastic system incoiporating 
feedback control forces and moments with time-delays 
in the state feedback is of an evident complexity, for a 
better understanding of the problem and of the present 
procedure, a simplified model has been adopted. 

DELAYED AEROELASTIC SYSTEM: 
STABILITY AND RESPONSE 

Some concepts related with the aeroelastic response and 
stability of the 1-DOF plunging airfoil in the presence 
of time delays between the sensing and the action of the 
actuator are presented next. 

A 1-DOF plunging airfoil is modeled as [1,4]: 

mh(t)+cli(t)+kh(t)=-L a (t) + L h (t)+ L c (t) (4) 
In the right hand side of this equation, the unsteady 
aerodynamic lift is represented as 

L a( t )= c iaP u M </>{t-a)^^-da + \pC La b 2 ii (5) 

The non-circulatory components of the unsteady 
aerodynamic load have been represented in terms of a 
convolution integral of the indicial Wagner’s function 
</>( r) , where the added mass is associated with the term 

jpC La b 2 h ■ 

In order to be able to highlight the implications of the 
nonlinearity in the delayed control, the structural and 
aerodynamic models are considered linear. In principle, 
with the exception of the unavoidable computational 
expense, the inclusion of these nonlinearities does not 
constitute a problem. For an aeroelastic model in which 
structural nonlinearities were included, see e.g. [4]. In 
addition, in Eq. [4] L h (r) denotes the external time- 
dependent load acting on the rigid wing counterpart and 
L c (r) denotes the nonlinear feedback control force: 

L c(t)= gp h (t - h ) ■ + ~h)~ g n y (t-t 3 ) (6) 

In the present work, the proportional (PFC) and velocity 
(VFC) feedback controls have been supplemented by a 
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nonlinear proportional feedback control (with delays 

*,•;* = 1,3). 

Denoting T = tU x /b and £ = h/b , the governing 
equation of the system with the nonlinear actuator 
control force with delay can be written in dimensionless 
form as 

r+2^f£'+(f )V 

v v J 

Structural terms 

v V— ' 

Aerodynamic terms 

+ g p (fff(r-r 1 ) + 2Qg v ff'(T-T 2 ) 

v v J 

Linear control terms 

-g„ c (t) 2 ^(r-r 3 )+ L h {r) 

v — > v ' v v- — ' 

Nonlinear control term External load 

The following dimensionless parameters have been used 
ju = m/ xpb 2 ; V =U ao /a> 0 b-,a> =® A /® 0 \a>l =k/m ; 
=c/2mw h ;g,, =g v /c,g p =g p /k-g nc = gnc b 2 /k . 

Evaluation of High-Order Aeroelastic Kernels 

Paralleling the procedure presented in [4], assuming a 
periodic external excitation of the form: 

/,,b)-X: (8) 

the high order kernels of the aeroelastic system can be 
derived. 

The identification of the n order aeroelastic kernels is 
based on a general input in the form given by Eq. (8) 
and on the extraction, for the generic term of w-th 

order, of the coefficients of n: -\ x< e • 

This procedure was detailed in [4], where the 
expressions of the first three Volterra kernels of 2-D 
lifting surfaces have been explicitly derived. 

As a remark, assuming a solution of the plunging 
displacement in the form %(t) = H lQL {s)X 1 e ST + ... , the 
first order kernel, // l0| (.v ) characterizing the open-loop 
1-DOF aeroelastic system can be represented as: 

H lOL {s)^ 2 + 2^fs + (fJ + j^{s) + \YY (9) 

In the present case we assume that the feedback control 
is represented by p(s) ; f)[s) can be one of the PFC, 
VFC, feedback control gains, or combinations of these 
(CFC), see [1], in conjunction with a nonlinear 
proportional control gain. 

Based on these assumptions, the linear first Volterra 
kernel of the closed loop system H lCL (s) is given by: 

H 1cl {s )=. tflOL (*)/[! + ^ 1 0I »/^)] (10) 

and is expressed explicitly as 


tfio»= { 2 + 2^ fs + (fj +^Ms)+it 2 

l-i (11) 

+ 9p{fJ e -*+2£s9 v fse-*] 

where the feedback gains are taken in absolute value. 
Usually, these gains are negative in the LQG/LQR 
design methodologies. 

For the specific case of the 1-DOF airfoil with linear 
control p{s) = (f )V sr > + 2Qg v fse~ ST2 . 

From (10), when p{s)= 0 , i.e. in the case of the open- 
loop, consistently H lQL = H lCL . 

From a mathematical point of view, the closed-loop 
nonlinear aeroservoelastic system, characterized by the 
first 7/| 0I (.S') and third H iOL (s l ,s 2 ,s 2 ) order kernels 

and by the feedback gain p{s) , can be seen as an open 
loop system described by the closed loop aeroelastic 
kernels // 1CL (s'), H 3CL (s l ,s 2 ,s 2 ) that are related to 

the kernels of the open loop system and its control 
gains. It should be mentioned that the second order 
kernel of the actual system is zero by virtue of the fact 
that no quadratic terms are involved in this system. To 
obtain the third-order closed loop nonlinear 

aeroservoelastic kernel we assume that 

the input can be expressed as L b (r) = . The 

output i §(t) can be written as: 

Z * Z>t M x V v 

u3) 

+ 3H 3cL {s 1 ,s 1 ,S}}X 1 X 2 X^ s,+S2 * Slir + others 

Substituting the expression of ^(r) as given by Eq. (13) 
in Eq. (7) and extracting the e s,T term and using 

the expression of (s t ) , the closed-loop third-order 
aeroelastic kernel is obtained 

^3Cl(' S ' 1’' S '2>' S '3) = 2g„c(f-) ^lCL^l^lCLfe^lCL^)^^ 

xH lCL {s l+ s 2 +s 3 )e-^ +S2+s ^ 

It is a general property of systems that all higher-order 
kernels can be expressed in terms of the lower-order 
kernels. In the absence of the quadratic term, see Eq. 
(7), the second-order kernel vanishes. Therefore, only 
the first and third order kernels have to be considered 
toward the determination of the nonlinear aeroelastic 
response and of the stability boundary. 

Stability Analysis 

For stability purposes, the aeroelastic system in the 
absence of external excitation, L b (r)= 0 , is considered. 
Without aerodynamic terms (that include time-lags), 
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and in the absence of control (i.e. g p = g v = 0), the 

system is dissipative with two finite stable characteristic 
roots (poles) on the left half of the complex plane. 
However, for the aeroelastic system with feedback 
delayed control ( t >0; j = 1,2,3 ) the two finite stable 

roots are supplemented by other finite stable roots 
(whose number depends on the aerodynamic model), 
and by an additional infinite number of roots due to the 
presence of e ST into the characteristic equation. 

The conditions that guarantee the stability of the 
delayed system, have been studied by Pontryagin [10], 
and applied toward the stability of time delayed 
feedback control systems by several authors. In the 
present aeroelastic analysis, Pontryagin’ s approach [10] 
in conjunction with the Stepan’s theorems [13] have 
been adopted. 

As proved in [12], the stability of delayed aeroelastic 
systems analyzed by using the concept of Retarded 
Functional Differential Equation (RFDE) depends on 
the presence of zeros with positive real parts of the 
characteristic equation, i.e. the presence of the p-zeros. 
For the stability evaluation, Eq. (11) can be written in 
characteristic equation form as: 

D(. S )=1/// Icl (.S')= 0, (15a) 

As a remark, in the absence of time delay, the following 
relation is valid 

D(s) = s 2 + 2Q f S + (f) 2 + J [o(,)+ 
+g p {f) 2 + 2 ^g v fs = o 

Since 1 + ^-[0(5) + ^-]> 0 , the stability conditions are 


obtained by imposing g y > -1 and g^ > -1 . 

Note that the characteristic roots (i.e. poles), of Eq. 
(15a) are of the form s - a + ico .Asa particular case, if 
g /; = g v = g , and r = Tj, the following relation holds 


valid: 


9 = 


1 + ^ 


M+i]fc)^ 2 +2^fe> + H 


! + 2^- 


(16a) 


It is possible to observe that for the uncontrolled system, 
g /; =g y = 0, the characteristic equation (15) has four 


finite stable poles in the complex plane that are obtained 
by solving the equation 


[l + A[q,(,) + i]j(E)V +2 ^(E) s + i| 


= 0 (16b) 


and all remaining poles are at a = -00 . For g ; , = g, = 00 , 


there are poles at s = 0;—l;— 82.36;— 12.89 and the 

remaining poles are at a = +00 . 

For equal time delays, via time transformation with 
respect to the delay, i.e. replacing sr => s , Eq. (15) can 
be further simplified as 


£$) = S 2 + 2C S fsr + (ffr 2 +^{S/t)S 2 

' / \2 - . (17) 

+ ys 2 +v p (fj? 2 e s +g v 2^f sre s = 0 

The stability of Eq. (17) will be studied via Stepan’s 
analytical method [8,13]. Following this method, upon 
denoting p x >...> p r >0 and CTj > . . . > a s = 0 , the 
non-negative real zeros of 

/?(©) = Re Z)(/®) = (-l)“ffl" +o(®") (18a) 

and 

S(co)=\mD(ico)=o{co n ), (18b) 

the trivial solution x = 0 of the system is exponentially 
asymptotically stable, if and only if 7?(o)>0 ; n = 2m 
( n is the order of the system and m is integer); 
S{p k ) ■*- 0 for k = 1, . . . , r and 


ZLi ^ s 8 n s (pk)= (“ O '" m ( ■ 1 ' 9 ) 

Similar conditions of stability are defined for systems 
where n = 2m + \, see [8]. For the present case, 
replacing o(i)5 2 => C(s)st and s=>ico, where c(k) 
(= F’(k)+/G(k)) is the Theodorsen’s function, and 
considering the real and imaginaiy parts of Eq. (17), we 
obtain: 


R(ct») = -o 2 +(y) --^-®?G(A:/r)to--^® 2 

+ g p (y f r 2 cos CO + 24 £ fg v rco sin CO 
5(®)= 2£ c yT co + a>rF(k /?) 

~9 P {ff r 2 sin ® + 2 ^ fg v rco cos co 


(20a) 


(20b) 


The trivial solution of Eq. (7) is exponentially 
asymptotically stable [8,23], if, and only if: 

9 P >-l (21a) 


CD - 


v. 9v 

9p 


< fiwf if - iff f + ifF OSfT + 2^ f Sine 
+ ^-(Gcoscr + Fsincr)} 


(21b) 


(21c) 


Herein, cr is the smallest positive zero of the equation 
* s ( cr ) = 2 ^fra + ^arF 

2 _ (22) 
-9p[f) T sm cr + 2^ f-g v rcr cos cr = 0 

where a e (0, y ) . The proof of Eq. (17) is given next. 
The inequality g > -1 is obtained from the condition: 

R(0)=(f) 2 f 2 +g/) (f) 2 ? 2 >0 (23) 

Considering the smallest positive root cr of Eq. (20b), 
one obtains that, ) > 0. co e (0, cr) if and only if 
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2 Qf® + jcoF-g p (ffrsinco 
+ 2Qy g v a> cos co > 0 

that yields 

g p {ff t sin co < 2£ ? f co + A coF + 2Q f g „ no cos co (25 ) 

Since the first two terms of Eq. (25) are always positive, 
it is possible to conclude that the condition 

r <24% Ifjp is required. In addition, for the smallest 

positive root cr , we can state: 

R{a) = -a 2 +(ffr 2 -jjarGa-jcr 2 

+ 9 P (f)^ 2 cos o- + 2£ % fg v rc sin cr < 0 
After some algebra manipulations 


(26a) 


g cos <j < 


V -1 + 2-Wu) 2 G _J_(.E^) 

\(D T ) U T \0) / LI \CO T f 


- 2 ^Mf)3r sinCT ] COSC 


(26b) 


and using standard trigonometric relationships, Eq. 
(26b) in conjunction with Eq. (25), rewritten in the form 


(if F 

+ 2^(ff)g vC rcos C r]s in ( T 


(27) 


reduces to Eq. (21c). The condition (19), where m = 1, 
is also fulfilled. 


time-delayed feedback control L c into Taylor series and 
omitting the terms containing the nonlinear function of 
the time-delay yields: 

L c {t)=-g v rh{t)+ (g„ - g p r)h{t)+g p h{t) 

+ g nc h 3 (t)-3g nc 7h 2 {t)h(t) + 0 (f 2 ) 

As a result, the aeroelastic governing equation is 
rewritten in the form 

mh(t) + ch(t) + kh(t) + g v rh(t ) - (g v - g p r )h(l) - g p h(t) 

~gj' 3 (/) + 3 g„ c rh 2 (t)h(t) = -C La pbU m j </>(t - a)h(t)da ( 2 
-\pC La b 2 h(t) + L c (t) 

and in dimensionless form as 

r + f? + (f ft = ~iiy(r - at'do ~ ^ 

-2^Qvf + kfQvf - + 9 p(f) 2 { W 

+ 9,,c{f?Z 3 - 3 9nc{f?rf~Z' 

Usually the gains are negative in the LQR design, and 
so are taken in absolute value. This yields the following 
approximated nonlinear equation: 

(l - 2^g„ f t-$' + y f + 2^g v f - 9p {f f rj? 

+ (f F (l + 9 P k + 9nc {ff fe 3 - 3?£ Y) W 


+ - 


f <^(r-af"da = 0 

J— 00 


The Aeroelastic Stability Chart 

The present approach for the determination of the 
stability domain of a delayed aeroelastic system has 
some analogies with Theodorsen’s method used for 
finding the flutter speed by plotting the real and 
imaginaiy parts of the flutter determinant in conjunction 
with consideration of a real co . The former approach 
reduces to the latter one in the case of zero time delays. 
The stability chart of the aeroelastic system described in 
Eq. (7) with respect to the feedback gains and the time- 
delay can be constructed using Stepan’s theorem and the 
D-subdivision method [12]. From the preliminary 
findings it appears that the time-delays play an 
important role. It is noted that the stability boundary 
depends dramatically on the velocity feedback control, 
especially in the case of the time delay. In addition, this 
implies that, in the presence of the delays, a small 
variation in the velocity feedback gain can expel the 
system from a stable to an unstable domain. 

HOPF-BIFURCATION ANALYSIS 


Preliminary findings related to the Flopf-bifurcation 
analysis of nonlinear feedback time delayed closed-loop 
aeroelastic systems will be presented next. 

Based on [8,24], assuming for sake of simplicity equal 
time-delays = r 2 = r 3 = r , expanding the nonlinear 


Stability examination 

The zero solution of Eq. (31), for the system in vacuum, 
is exponentially asymptotically stable [8,24] if 

9 P > -1 (32a) 

l-2^g„t?>0; (32b) 

2^f + 2^g v f-g p (f) 2 f>0 (32c) 

The FIopf-Bilurcation (HB) is present at 


9; 




; g,. <- 


g„ >-l (33a, b,c) 


(l + g v ) fr- 2C f - r . 

and it is supercritical or subcritical if g (!C <0 or 
g„ c > 0 , respectively. The proof, paralleling that 
presented in [8] is given next. 

The nonlinear Eq. (31) can be linearized about g = 0 

(1 - 24 $ v f r)r + y f + 24 & f-g P (t) 2 y 

+ (t) 2 (i + g^ = o 

and its characteristic polynomial is given by 

(l - 24 # v f F> 2 + [ 2 ^ f + 2Qg v f - Qp [ffz\ 


(34) 


(t) 2 ( 1 + 9p)=0 


(35) 


Assuming that g v < and g v = ~^9 P -1 , then, 

as a necessary requirement for the Hopf-bifurcation, we 
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have the condition of a pair of pure imaginary roots of 
the characteristic equation 


s l2 =+ij3 where j3 = 


1 + 9 , 




= r -r 2 g 

m I- 


(36a, b) 


If V = 0 then R = ■ 

i~r 2 g p 

Upon defining the new variables y 1 = £ and 
y 2 = £,’ / P the Poincare normal form is obtained: 


y I 


' 0 -p 

)’i 


0 

y\_ 


_-P 0 . 

_y 2 _ 

+ 

./30T1 3 + f 2 \y\y 2 _ 


where 

/3o = -^r(i+g P )g n c ( 38a ) 

/2i = f(i+g P kc=-3^/ 30 (38b) 

The type of HB occurring at the critical time-delay can 
be determined by using the center manifold theorem. 
Specifically, the sign of the quantity L , where 

L = i/21 (39) 

determined at the critical time -delay, defines the benign 
flutter boundary (supercritical HB) and the catastrophic 
one (subcritical HB), for the cases L < 0 and L > 0 , 
respectively. This implies that 

/ 21 < 0 => g„ c < 0 , supercritical HB (40a) 

/ 21 > 0 => g nc > 0 , subcritical HB (40b) 

g„ c = 0 degenerated HB. (40c) 


RESULTS AND DISCUSSION 


The parameters for the simulation are presented in Table 1. 


Table 1. Airfoils and flow parameters 


l-DOF Plunging Airfoil 

6 = 1ft 

p x = 0.03 18 slugs/ ft J 

^ = 0.008 

c = 2 mco£ 

m = 10 slugs/ft 

k = co 2 m 

a> h = 60 rad/s 

C La =2K 


In Fig. 1 the effect of the time delay on the first order 
kernel is presented for selected values of the 
proportional and velocity feedback gains and for 
selected values of the dimensionless frequency , 

respectively. In Figs. 2 and 3, 3-D and contour plots of 
the first and third order kernels with the variation of the 
frequency and of the time delay are depicted for two 
values of the proportional and velocity feedback gains 
and two values of the time delay, respectively. 

The stability chart of the aeroelastic system described 
by Eq. (7) with respect to the feedback gains and the 


time-delay is constructed using the second Stepan’s 
second theorem and the D-subdivision method. 


0.015 

A g p =g v =io 4 

/ 1 — no delay 

0.0125 

/ 

^ 0.01 

3^ 


J 0.0075 

V\\ \\ Time delay 

v\\ A T- 0:0.01; 0.02; 0.03; 0.04; 0.05; 0.06 

0.005 

yC\\ 

0.0025 


0 




0 20 40 60 80 100 

Frequency, to 



Time delay, x 


Fig. 1 First order aeroelastic kernel. Effect of the time 
delay, a) Hi vs. a> . Influence of the time delay; b) /// vs. 
t. Influence of the frequency. 

The method of D-subdivision is applied for determining 
the condition under which the quasi-polynomial D c (s ) 
has no p-zeros. 

As remarked in Ref. [3], since the quasi-polynomial is a 
continuous function of its parameters we can construct 
the subdivision of the coefficient’s space by hyper- 
surfaces, the points of which are quasi-polynomials with 
at least one imaginary root. In addition, as it has been 
proven in [12], with the variation of the quasi- 
polynomial parameters the number of p-zeros may 
change only by passage of some zeros through an 
imaginary axis, and that for all points of every domain 
of D-subdivision the number of quasi-polynomial p- 
zeros will be the same. 

The region in the {g , g v } parameter space where the 

roots of the characteristic equation of the system have 
zero real part has been represented in Fig. 4 in which the 
values of the geometrical parameters have been locked. 
Other stability regions can be drawn in the {g^,g v j 

parameter space, but these are not of practical 
importance. As clearly presented, the time-delays play 
an important role. It is noted that the range of stability 
for g is much larger than that of g r , implying that the 

stability boundary depends dramatically on the velocity 
feedback gain, especially in the case of the time delay. 
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g = g = 

-‘p 


1.0 - 10 4 



0 0.2 0.4 0.6 0.8 1 

Time delay, x 


Fig. 2 3-D and contour plots of the first order kernel. Variation of the frequency vs. time delay for two values of the 

proportional and velocity feedback gains. 
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Time delay, x = 0.1 



-40 -20 0 20 40 

032 Time delay, x = 0.25 

g =g =2.5 -10 4 

~ J p V 



-40 -20 0 20 40 


®2 


Fig. 3 3-D and contour plots of the third order kernel. Variation of the frequency vs. the time delay for two values of 

the time delay. 
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In addition this implies that, in the presence of the 
delays, a small variation in the velocity feedback gain 
can expel the system from the stable domain to the 
unstable domain. On the other hand, the stable 
parameter space is the complete positive quadrant of 
jg ,g v } parameter plane if no time delays are present, 

independent of the values of the feedback gains. 


■nf 



Fig. 4 Stability charts of the full and approximate 
(Taylor’s expansion) closed-loop delayed nonlinear 
aeroelastic system. 

Whereas the linear analyses are able to predict the 
stability boundary, the nonlinear analyses provide an 
insight on the character of the stability boundary. To 
this purpose an auxiliary stability chart (see Fig. 4), 
corresponding to the statement of Eqs. (33), enables the 
prediction of the type of FIopf-Bifurcation (i.e. 
supercritical or subcritical HB). This auxiliary plot for 
small time-delay, supplies full information related to the 
FIB but, due to the simplicity of the Taylor’s series 
expansion used, provides other types of instability. For 
this reason the auxiliary plot is used only to determine 
the FIB conditions and not the analytical examination of 
the full nonlinear system. 


■id 5 



Fig. 5 3-D stability chart of the plunging airfoil. 
Effect of the time delay and of the control gains. 


On the basis of the FIB condition and of Fig. 4, it can be 
seen that the FIB occurs on the 3-D surface of the 
parametric domain depicted in Fig. 5. Therefore, if the 
parameter point is situated in the inner domain, the 
system is stable, and if the parameter point is in the 
outer domain, the system reaches a stable or unstable 
LCO if the nonlinear control gain is negative or 
positive, respectively. In Fig. 6 the nonlinear closed- 
loop aeroelastic responses of the 1-DOF airfoil for 
g„ c =-10 (implying a stable LCO, i.e a supercritical 
FIB) and g /!C = 10 (implying an unstable LCO, i.e. a 
subcritical FIB) are presented. These results, obtained 
via both Volterra series and numerical integration, are in 
agreement with the analytical predictions of Eqs. 
(40a, b). 


stable LCO 



4 unstable LCO 



Fig. 6 Stable and unstable LCO for the plunging airfoil 
with active feedback delayed control law (r = 0.0 1) . 

CONCLUSIONS 

This paper constitutes a basic step towards a better 
understanding of a number of issues related to the 
nonlinear aeroservoelasticity of 2-D and 3-D aircraft 
wings and their nonlinear feedback control featuring 
time-delays. The model presented here, based on the 
determination of higher order Volterra kernels, can be 
used toward the determination of the nonlinear 
aeroelastic response and flutter boundary of 
open/closed-loop aeroelastic systems. The effects of the 
control parameters and of the delay in the state feedback 
on the flutter boundary and on the stable/unstable LCO 
have been highlighted. The importance of the time delay 
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on the nonlinear feedback control has been emphasized. 
With the incorporation of structural and aerodynamic 
nonlinearities in the system, or even in the feedback 
loop, more complex phenomena (LCO and chaos) are 
expected to occur. 
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